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Red star-forming and blue passive galaxies in clusters 
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ABSTRACT 

We explore the relation between colour (measured from photometry) and specific star for- 
mation rate (derived from optical spectra obtained by SDSS DR4) of over 6000 galaxies 
{Mr ^ —20.5) in and around (< 3 r2oo) low redshift (z < 0.12) Abell clusters. Even though, 
as expected, most red sequence galaxies have little or no ongoing star formation, and most 
blue galaxies are currently forming stars, there are significant populations of red star-forming 
and blue passive galaxies. This paper examines various properties of galaxies belonging to 
the latter two categories, to understand why they deviate from the norm. These properties 
include morphological parameters, internal extinction, spectral features such as EW(H5) and 
the 4000 A break, and metallicity. Our analysis shows that the blue passive galaxies have 
properties very similar to their star-forming counterparts, except that their large range in 
equivalent width indicates recent truncation of star formation. The red star-forming galaxies 
fall into two broad categories, one of them being massive galaxies in cluster cores dominated 
by an old stellar population, but with evidence of current star formation in the core (possibly 
linked with AGN). For the remaining red star-forming galaxies it is evident from spectral in- 
dices, stellar and gas-phase metallicities and mean stellar ages that their colours result from 
the predominance of a metal-rich stellar population. Only half of the red star-forming galaxies 
have extinction values consistent with a significant presence of dust. The implication of the 
properties of these star-forming galaxies on environmental studies, like that of the Butcher- 
Oemler effect, is discussed. 

Key words: Galaxies: clusters: general; galaxies: evolution; galaxies: fundamental parame- 
ters; galaxies: starburst; galaxies: stellar content. 



^ ' 1 INTRODUCTION 



If, as is widely believed, large-scale structures form in the Universe 
by means of the hierarchical assembly of dark haloes, one would 
expect the most massive galaxies to be the youngest, as they would 
be formed from mergers over time. However, the brightest galax- 
ies, elliptical in morphology, typically found at the centres of rich 
clusters, appear to be dominated by old stellar populations. Galax- 
ies with current star formation, on the other hand, have been shown 
to be disc-dominated, mostly occurring in regions of lower density. 
Such relations between the properties of galaxies and their envi- 
ronment were first pointed out using morphology to characteris e 
galaxies l lOemlei|[T974l : iMehiick & Sargent|[l97l , lDressleJI 19801) , 
but, in later studies, other photometric properties such as colour 



(e.g.lKauffmann et alj20"033 : lBalogh et alj2004al : IWeinmann et all 



l200alCooper et al.ll2007h " have been shown to correlate with envi- 
ronment as well. 

In addition to understanding the formation of dark matter 
haloes, the study of galaxy formation and evolution clearly re- 
quires detailed knowledge of the various physical processes that 
affect baryons, such as mergers, tidal interactions, stripping, cool- 
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ing and feedback. Observations of the relevant galaxy features and 
their dependence on environment are very important in our com- 
prehension of the relative importance of these processes, a sub- 
ject that remains a n issue of considerable disagreement (see, e.g., 
IOemleretalj2009l) . 

The availability of large-scale spectroscopic surveys, 
such as the Sloan Digital Sky survey (henceforth SDSS), 
has resulted in structural parameters of galaxies being 
largely replaced by spectroscopic properties, in the explo- 
ration of the impact of th e global an d local environment 
on the n ature of a gal a xv (jSalogh et alj l2004bl: iHaines et al.l 
20063,"; Pimbblet et al. 2006': 'Porter & Raychaudhury' l200a 
Weinmann et al, 2006; Porter & Raychaudhury 2007; Porter et al.l 

in preparation! , among 



20081 : iMahaian. Raychaudhury & Pimbbled 



others). Compared to broad-band colours, spectroscopic param- 
eters such as integrated star formation rate (SFR), specific star 
formation rate (SSFR or SFR/M*), equivalent width (EW) 
and the 4000 A break (D,i4000), in principle, give more detailed 
information about the underlying stellar populations in galaxies, 
but one cannot expect these parameters to behave in the same way 
with environment as the global photometric properties. 

Various factors, including the presence of dust, and the well- 
known degeneracy between the effect of age and metallicity, pre- 
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vent a straight-forward correspondence between the s tar formation 
history of a g alaxy and its broad-band colours (e.g., iFab^l 19731 : 
ICalzettilll997h . As a result, the use of broad-band colours and of 
spectroscopic indicators might yield different classification for a 
galaxy on an evolutionary sequence. In order to compare obser- 
vational studies where different photometric or spectroscopic pa- 
rameters have been used to define sub-samples, one would like to 
ask, for instance, for what fraction of galaxies the correspondence 
between morphology, colour and star formation properties fails to 
hold, and whether such galaxies could be identified from other ob- 
servable properties. Are these unusual galaxies more likely to occur 
in certain kinds of environment? 

In order to address some of these questions, we choose, for 
this study, galaxies in and around rich clusters in the low redshift 
Universe. Galaxy clusters are well suited for such an analysis, since 
from their cores to their infall regions, they provide a wide range of 
environments, with varying space density and velocity dispersion. 
The aim of this work is to examine how well galaxy properties 
obtained from two mutually independent techniques, namely pho- 
tometry and spectroscopy, correlate with each other, and to explore 
the possible causes of discrepancy. 

The outline of this paper is as follows: in the next section we 
define the galaxy sample, and the galaxy properties to be used here. 
We plot the specific star formation rate of these galaxies against 
colour, and identify the four classes of galaxies that we will deal 
with: red sequence, red star-forming, blue star-forming and blue 
passive. In ^ we examine the various galaxy properties, derived 
from photometry and spectra, that we use in this analysis. The in- 
sight these properties provide into the nature of blue passive and 
red star-forming galaxies are discussed in [j?] and the main con- 
clusions are summarised in Sj5] Throughout this work we adopt a 
cosmology with flm = 1 & = and h — 70 km s^^ Mpc^^ 
for calculating distances and absolute magnitudes. We note that in 
the redshift range chosen for this work (0.02 ^ 2 ^ 0.12), the 
results are insensitive to the choice of cosmology. 



2 THE OBSERVATIONAL DATA 

We u se the spectroscopic galaxy da ta provided by the SDSS 
DR4 jAdehnan-McCarthv et alj[200 ?). The SDSS consists of an 
imaging survey of tt steradians, mainly in the northern sky, in five 
passbands u,g,r,i and z. The imaging is done in drift-scan mode, 
and the data is pro cessed using the photometric pipeline PHOTO 
jLuDton et al .l200lh specially written for SDSS. The spectra are ob- 
tained using two fibre-fed double spectrographs, covering a wave- 
length range of 3800-9200 A. The resolution A A/A varies between 
1850 and 2200. 

We use the galaxy magnitudes and the corresponding galac- 
tic extinction values in the g & r SDSS passbands from the 
New York Universi ty Value added galaxy Catalogue (NYU-VAGC; 
iBlanton et al.l200^ . We fc-co rrect these magnit udes to z = 0. 1 (me- 
dian redshift of the sample, lYang et ^|2007|) . We only consider 
galaxies with absolute magnitude ^ — 20.5 (this corresponds to 
the apparent magnitude limit for the SDSS spectroscopic catalogue 
at the median redshift of our sample, z — Q.l). 



2.1 Cluster membership 

We take all Abell clusters dAbell et al.|[l989l) , in the redshift range 
0.02< z ^0.12, that lie sufficiently away from the survey bound- 
ary so as to be sampled out to 6 h^g Mpc from the cluster centre. 



We constrain our sample to moderately rich clusters by considering 
only those that have at least 20 galaxies (brighter than the magni- 
tude limit, Mr ^ -20.5) with SDSS spectra within 3 h:^^ Mpc of 
the cluster centre. This leaves us with a sample of 119 clusters. 

We aim to study the effect of present and recent star for- 
mation on the observed properties of galaxies over a wide range 
of environment, and so we seek to examine galaxies not just in 
the cores of the clusters, but also on the outskirts. Detailed stud- 
ies of the velocity structure in the infall regions of clusters show 
that the turnarou nd radius of clusters extends out to almost twice 
the virial radius jPiaferio & Gelleij|l997l : iRines & Diaferioll2006l: 
iDtinner et alj|2007[) , so we chose to examine radial trends of the 
properties of galaxies out to a projected radius of i? = 6 h^^ Mpc 
of each cluster centre. Within this radius, we pick galaxies within 
the following velocity ranges about the mean redshift of the cluster: 



3000 km s"^ 
2100 km s"^ 
1500 km s"^ 



if 7? «; 1 Mpc 

if 1 < i? ^ 3 Mpc 

if 3 < i? ^ 6 Mpc 



(1) 



Adopting these differential velocity slices helps in reducing inter- 
lopers, in the absence of a rigorous procedure for assigning cluster 
membership, and hence is well suited for a statistical study like this. 
Note that other popular methods, such as the 3-(t„ velocity disper- 
sion cut, are prone to include galaxies at large projected radius from 
the cluster centre, and our velocity cuts are consistent w ith the caus- 
tics in velocity sp ace used by the CAIRNS survey (e.g. lRines et al.l 
I2OO3I: iRines&Diaf eriQi2006i) . 



2.2 Cluster parameters 

We find the velocity dispersion (a^) of each cluster, from the 
iden tified memb ers, using the bi-weight statistics (ROSTAT, 
iBeers. Flvnn & Gebhardt 1990) . The characteristic size of the clus- 
ter (the radius at which the mean interior over-density in a sphere of 
radius r is 200 times the critica l density of the Universe) i s calcu - 
lated using the relation given bv lCarlberg. Yee&Ellingsoril ( ll997h : 
7'200 = \/3cr„/10-ff (2). Our final sample comprises of 119 Abell 
clusters (z < 0.12) in which there are 6,010 galaxies within 3 r2oo 
of the cluster centres. 



2.3 Galaxy parameters 

In this work, we u se various quantiti e s de rived from the SDSS 
optical spectra by iBrinchmann et alj ( 12004) (B04 henceforth), 
based on SDSS DR4 c atalogujE From this and related work 
jKauffmann et al. Il2003lj), we use specific spectral indices (emis- 
sion corrected equivalent width and [Nelll] contamination cor- 
rected D„4000). 

Two of the galaxy parameters we use to begin with, namely 
the specific star formation rate (SSFR, or SFR/M*) and stellar 
mas s (M*), are derived using the SDSS galaxy spectra by B04 
and iKauffmann e t al. ( 2003b) respectiv ely. B04 construct a grid 
of ~ 2 X lO'^ spectral syn thesis models dSruzual & Charlot|[T993l: 
ICharlot & Longhettill200ll) over four relevant parameters spanning 
a range in metallicity (Z) of —1 < log Z/Zq < 0.6, ionization 
parameter (U) —4.0 < log U < —2.0, total dust attenuation (rv) 
0.01 <Tv < 4.0 and dust-to-metal-ratio (0 0.1 < C < 0.5. The 
attenuation by dust is treated using a physically motivated model 



^ http://www.mpa-garching.mpg.de/SDSS/DR4/Data/sfr_catalogue.html 
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Figure 1. The distribution of the gal axies used in this work in the plot of 
the two line-ratios, recommended by [Baldwin, Phillips & TerlevichI jl98lL 
(BPT)), that are often used to discrimi nate between star-forming (and com- 
posite) galaxies (below the red line; iKauffmann et al.l i2003cl) ) and the 
AGN-dominated galaxies (above). All galaxies in our sample, except the 
'unclassified' , are plotted here. 



of lCharlot&Fdl l l2000h . which provides a consistent model for 
ultra-violet (UV) to infra-red (IR) emission. 

B04 compare the modelled line-ratios for each model in their 
grid directly with those in the observed spectrum, and, for each 
galaxy, construct a log-likelihood function for each point on the 
model grid, summed over all lines observed in the spectrum of that 
galaxy. It is worth mentioning here that since they study relative 
line fluxes only, their analysis is insensitive to stellar ages, star for- 
mation history and the relative attenuation by dust in the molecular 
gas clouds and the interstellar medium. 

B04 derive their values for star formation rate (SFR) and spe- 
cific star formation rate (SFR/M*) using a Bayesian technique, in- 
volving the likelihood functions mentioned above, using the en- 
tire spectrum of each galaxy (it is worth noting that the spectrum 
comes from a fibre of diameter 3"). It is shown that their work yield 
better estimates of the SFR and SFR/M* than any single spectro- 
scopic SFR indicator, such as equivalent width (EW) of the Hc« or 
O II lines, which are used extensively in the literature. Furthermore, 
the combination of several line indices enables one to separate the 
galaxies with a significant AGN component from those with emis- 
sion entirely due to star formation. The derived parameters (SFR 
and SFR/M*) are then corrected for aperture biases using their pho- 
tometric colours. Out of the three statistical estimates (mean, me- 
dian and mode) for the probability distribution function (PDF) of 
SFR and SFR/M* derived for each galaxy, we use the median of the 
distribution in this work for each quantity, since it is independent 
of binning. 

Star-forming galaxies can be separated from the AGN- 
dominated ones on the basis of the so-called BPT diagram 
tealdwiiT Phillips & TerleviclJll98lh . namely a plot of the line ra- 
tio [OIII]/H^ against Hq/[NII], Fig. [T] for galaxies with all four 
emission lines. B04 divide all galaxies into four classes (star form- 
ing, AGN, composite and unclassified) using the BPT diagram. For 



those galaxies whose spectra do not show all the above mentioned 
emission lines, (called 'unclassified' galaxies), and those that are 
classified as AGN-dominated, B04 adopt an indirect method for es- 
timating the SFRs. They use the relationship between the D„4000 
and SFR derived for the high S/N (^ 3) star-forming galaxies 
(SFRe) to estimate the SFR of the AGNs and low S/N galaxies 
(SFRti; see section 4.1 of B04 for details). Using the entire spec- 
trum also makes it possible to obtain an independent estimate of 
stellar mass M*, which further constraints the estimated SFR/M*, 
making it a more reliable parameter for evaluating star formation 
properties of galaxies, compared to other luminosi ty based pa- 
rameters employed elsewhere (e .g. [Lewis et al. 2002; Gomez et al.l 
l2003h . lKauffmann et al.l ( l2003ij) adopt similar methodology as B04 
for estimating the stellar mass M* . Their models are characterised 
by star-to-light ratio, dust attenuation and fraction of mass formed 
in a major burst of star formation (the burst mass fraction). They 
then compute the total mass in stars by multiplying the dust- 
corrected and fc-corrected z-band luminosity of the galaxy, by 
the estimated z-band mass-to-light ratio of the best-fitting model, 
found using the Bayesian technique described above. 

Several authors have shown that the SFR of a galaxy scales 
with its stellar mass (e.g. Noeske et al. 2007). Since the reason for 
this correlation is yet unclear, and beyond the scope of this work, 
we employ the specific star formation rate SFR/M* as the more 
reliable measure of the star formation activity of galaxies. The val- 
ues of SFR (and SFR/M*) derived for the galaxies with emission 
lines, i.e. all but those categorised as 'unclassified', are reliable 
(see Fig. 14; B04), because the classification requires S/N> 3 in 
all the 4 lines involved in the BPT diagram (i.e.. Ha, H/j, [OIII] 
and [Nil]). But for spectra of galaxies without emission lines, or of 
those that are unclassified (most of which are passive), the uncer- 
tainties are relatively higher. We chose not to select galaxies based 
on the confidence interval values derived from the PDF of the SFR 
and SFR/M* parameters of individual galaxies, because excluding 
galaxies with higher uncertainties in SFR/M* (or SFR) distribution 
would preferentially exclude a significant number of passive galax- 
ies and thus bias our samples. 

2.4 Four categories of galaxies 

Galaxy properties measurable from photometry or spectra, such 
as colour, SFR or SSFR (SFR/M*) have been extensively used in 
studies of the ir depe ndence on the environm ent of the galaxy (e.g . 
'Balogh et al." 2004b; Haines et al. 2006a b; Pimbblet et al. 200^ 
iWeinmann et al..,2006; Porter & Ravchaudhurv. 2007; Porter et al 



I2OO8I ; Mahaian. Ravchaudhurv & Pimbblej in preparationh . or on 
other variables representing mass, such as the stellar mass 
of the galaxy |K auffmann et al. 2003b) or its circular velocity 
jGraves et al.l2007l) . 

In this work, we study the colour and SSFR as indicators of the 
evolutionary stage of a galaxy. In Fig.|2]we plot the colours of all 
galaxies found within 3 r2oo of the centres of all the clusters in our 
sample against their specific star forma tion rates. In a pion eerin 
study of evolution of galaxies in clusters. [Butcher & Oemlerl ( ll98' 
had adopted a colour cut to define blue galaxies ?& B — V — 0.2 
mag lower than the cluster's red sequence. Here we stack all the 
cluster galaxies together, and fit a red sequence to galaxies which 
has a slope of —0.034 with an absolute deviation of 0.080. Since 
this slope is negligible, we adopt a null slope red sequence at {g — 
r)"'^ _ I jjQjg jjjjf using the actual red sequence fit does not 
change our results at all. The blue galaxies are thus defined to have 
(5 ^ r)"'^ ^ 0.8. This classifies ~25% of all galaxies as blue. 
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Figure 2. The relation between specific star formation rate (SSFR or 
SFR/M* ) for all galaxies found within 3 r200 of the centre of each galaxy 
cluster in our sample. The contours represent 150, 100, 50 and 25 galax- 
ies respectively. The horizontal line marks the colour cut adopted by most 
Butcher-Oemler effect studies, i.e. 0.2 mag bluer than the red sequence: this 
labels ~ 25% of all galaxies as blue. To match this cut, we adopt a SSFR 
threshold, log SFR/M*= -10.16 yr^^ such that 25% of all the galaxies 
are classified as star-forming. We thus categoiise our galaxy ensemble into 
four classes: (clockwise from right) red star-forming, blue star-forming, blue 
passive and red sequence galaxies respectively. 



Table 1. Number of galaxies in the four categories 





red SF red sequence 


blue passive 


blue SF 


Clusters 


539 3998 


538 


935 


Entire SDSS 


24039 163904 


44877 


113909 



Correspondingly, we choose the SSFR threshold at log 
SFR/M* = -10.16 yr"\ defining the 25% of all galaxies with a 
liigher value of SSFR as star-forming. We note that our results do 
not qualitatively change even if the SSFR threshold is changed to 
the 50-percentile value (log SFR/M*= -10.45 yr"^). We thus clas- 
sify the galaxies in the 4 quadrants of Fig. |2] as red star-forming, 
red sequence, blue passive and blue star-forming, anti-clockwise 
from top right. Table[T]gives the number of galaxies in our cluster 
sample in each category, and also the corresponding numbers in the 
whole of SDSS DR4 spectroscopic catalogue in the same redshift 
and magnitude range as our cluster galaxy sample. 



2.5 Spatial distribution of galaxies in clusters 

Having divided the galaxies into four classes, we consider their ra- 
dial distribution within the clusters to examine the possible influ- 
ence of local environment. Fig.[3]shows the distribution of galaxies 
in all four categories as a function of distance from the correspond- 
ing cluster, scaled by r2oo- For each category, the curve has been 
normalized by the total number of galaxies in that category. 

As expected, the fraction of the passive red sequence galax- 
ies increase towards the cluster centre. The red star-forming galax- 
ies closely follow the distribution, but their fraction inside r2oo 
does not increase as steeply as that of the red sequence galax- 
ies, suggesting a different evolutionary path. Here it would have 
been very interesting to examine the detailed difference between 
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Figure 3. The radial distribution of galaxies in the four categories defined 
in il2.4l is plotted as a function of cluster-centric radius, scaled by r200- 
Each distribution is normalized to unity. The Poisson error on each point 
is smaller than the size of symbols. The fraction of red sequence galaxies 
(black), as expected, increases toward the cluster centre, and the red star- 
forming galaxies ( red) follow them closely. The blue star-forming galaxies 
(blue) show a decline within r200, while their passive counterparts (green) 
show a constant fraction out to 3 r2oo, except for a mild increase ~ 2r 200 ■ 



these two distributions as a function of halo mass, since in both 
clusters and groups the morphology and star formation properties 
seem to vary between dwarfs and giants (Khosroshahi et alJ I 20041 : 
Miles, Raychaudhury, & Russell 2006), and/or central and satellite 
galaxies dAnn, Park. & Choill2008l) . indicating different evolution- 
ary histories. But in this work, over the redshift range of our sam- 
ple, we have to be confined mostly to M« and brighter galaxies 
(Af^:^-20.5). 

At the other extreme, the blue star-forming galaxies show a 
constant distribution beyond r2oo, dropping steeply in the cluster 
core. This is the consequence of the observ ed lack of neutral gas 
in ga l axies in the cores of rich c lusters (e.g. iGiovanelli & HavnesI 
ll985l : lRavchaudhurv et al.lll997l) . 

However, the blue passive galaxies, which also have a more 
or less constant radial distribution over the entire range, do 
not show this rapid drop in the cores. This is consistent with 
the observation, discussed below in some detail, of nuclear star 
formation in a substantial fraction of br ightest cluster galaxies 
teildfell et al] I2OO8I ; iReichard et al] |2009|) . which will cause the 
innermost bin to be significantly higher for this distribution than 
that for the blue star-forming galaxies. It is interesting to note 
that the fraction of blue passive galaxies also show a minor in- 
crement around 2r2oo- Other studies have shown that in clus- 
ters, particularly those that are fed by filaments, there is an en- 
hancement of nuclear star formation in the infalling galaxies be- 

I 1' 1 I '•^-^ II — — 1 

yond r2oo (e.g. Porter & Raychaudhurv 2007; Porte r et al.M2008l : 

iMahaian. Ravchaudhurv & PimbbleLin preparatior;) . 



3 DO THE GALAXY PROPERTIES DERIVED FROM 
PHOTOMETRY AND SPECTROSCOPY AGREE WITH 
EACH OTHER? 

One of the objectives in this paper will be to examine the hypothe- 
sis that the galaxy populations can indeed be divided, on the basis 
of a single parameter, into two classes: passive and star-forming, 
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and that the incidence of red star-forming or blue passive galax- 
ies is merely due to the scatter in the measured parameters within 
reasonable errors. As we show below, this hypothesis fails to hold. 

3.1 The effect of extinction and abundance 

Star-forming regions in a galaxy are often surrounded by dust, 
which can significantly affect quantitative analyses based on opti- 
cal light. Conversely, the presence of dust in a galaxy would imply 
the presence of star-forming regions therein. In the absence of in- 
frared observations, which are often used to quantify the effect of 
dust in star-forming systems (e.g. Calzetti 1997), optical photome- 
try can also be used to constrain the internal extinction in a galaxy 
(also see Fig. |9}. 

In one such attempt, iKauffmann et alj ( l2003bh have estimated 
the degree of dust attenuation in the SDSS z-band, by comparing 
model colours to the measured colours to estimate the reddening for 
each galaxy. The shape of the attenuation curve from observations 
is found to resemble a power law with slope ■-^ — 0.7 over a wave- 
length range from 1250-8000 A jCalzetti et aT]|l994h . A standard 
attenuation curve, of the form t\ oc A""'^, is then extrapolated to 
obtai n the extinction in the SD SS z-band. We use the Az values 
from IKauffmann eTall ( l2003bi) . expressed in magnitudes, in this 
section. The typical la error on the estimated is ~ 0.12 mag. 
They also note that comparison of g—r and r—i colours yield sim- 
ilar extinction values for most galaxies if the colours are corrected 
for nebular emission. 

The reddening of a galaxy is the result of a combination of 
several factors, involving the presence of dust and/or metals in the 
interstellar medium and of old stars. In our attempt to interpret the 
apparent evidence of star formation in the fibre spectra of some red 
galaxies ('red SF'), and the low values of SFR in the spectra of 
some blue galaxies ('blue passive'), we plot, in Fig.|4] the distribu- 
tion of internal extinction Az of the galaxies in the four populations 
as classified in Fig.|2] 

As expected, the 'red sequence' galaxies, most of which are 
passively evolving massive galaxies in the cores of clusters, show 
very little extinction (negative values of Az being interpreted as 
zero here). However, even among passive galaxies, there is a small 
fraction that constitutes a tail extending to non-negligible extinction 
values. The blue star-forming galaxies, on the other hand, show a 
clear sign of attenuation going upto 1.9 mag in the z-band, with 
most galaxies having extinction values around 0.6 mag for the blue 
star-forming galaxies and 0.35 mag for the blue passive galaxies 
respectively. 

In contrast, the red star-forming galaxies show a mix of at least 
two different populations: about half of them are unattenuated (sim- 
ilar to the red sequence galaxies), and the other half contributing to 
the skewed high extinction end of the distribu tion. We employ the 
KMM algorithm jAshman. Bird & Zepj[l994l) to quantify the ex- 
istence of bimodality in our data. The KMM algorithm fits a user 
specified number of Gaussian distributions to the dataset, calculates 
the maximum likelihood estimate of their means and variances, and 
assesses the improvement of the fit over that provided by a single 
Gaussian. The KMM likelihood ratio test statistics (LRTS) used 
here is a measure of improvement in using two Gaussian distribu- 
tions (since we are testing bimodality) over a 1-mode fit. 

For the distribution of extinction, in the red star-forming 
galaxies the KMM test gives the value of likelihood ratio (LRTS) 
to be 75, when using a bimodal fit to the data than using a unimodal 
fit, with nil probability of the null hypothesis being satisfied. The 2- 
modes of the internal extinction distribution are found to be centred 
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Figure 4. The distribution of internal extinction among our cluster galaxies 
in the four quadrants as specified in Fig. |2] Each histogram is normalized 
to unity. It is clear that in general more 'blue' galaxies are prone to internal 
reddening than the 'red' ones, with star-forming (blue dotted) and passive 
blue galaxies (green dot-dashed) showing slight offsets. Interestingly, while 
< 45% of the red sequence galaxies {black solid) show signs of some ex- 
tinction, the red star-forming galaxies {red dashed) have a bimodal distri- 
bution in extinction values (see text). The grey lines represent the same dis- 
tributions as above, but for > 300, 000 galaxies found in the entire SDSS 
DR4, divided in the four quadrants in the same way as described for the 
sample of cluster galaxies studied here. Here, negative values of Az indi- 
cate the absence of extinction. 

around Az=-0.06 (read zero extinction) and 0.72 mag respectively 
for the sample of cluster galaxies. Interestingly, we find that even 
by changing the selection criteria for defining the red star-forming 
galaxies to a higher value of log SFR/M* = -10 yr^^, the KMM 
test still yields statistically significant result in favour of a bimodal 
fit to the data (LRTS=21 at a significance level 5 x lO"'^). 

Along with the distribution of extinction values for the clus- 
ter galaxies, we also show the corresponding distributions for 
>300,000 galaxies, with the same magnitude and redshift range 
as our cluster sample, drawn from the entire SDSS DR4 spectro- 
scopic catalogue (grey thin lines). For this much larger sample, the 
histograms of extinction values are not significantly different from 
those for the cluster sample, except for the red star-forming galax- 
ies, where cluster galaxies seem to have relatively fewer galaxies 
with high extinction values than similar galaxies in the larger all- 
SDSS sample. 

3.2 H« and the 4000A break 

In order to probe the presence of young (< 1 Gyr old) stellar pop- 
ulations among galaxies with an unusual combination of SFR/M* 
and colour, we plot the distributions of the EW of the absorp- 
tion line, and of the D„4000 values in Figs. |5] and |6] respectively. 
The D„4000 is the strongest discontinuity occurring in the optical 
spectrum of a galaxy, mainly due to the presence of ionised met- 
als. On the other hand, absorption is detected after most of the 
massive hot stars have finished evolving on the main sequence, i.e. 
at least 0.1-1 Gyr after a starburst is truncated. Thus the Hf EW is 
a measure of the age of the youngest stellar population in a galaxy, 
and has been extensively used to e stimate the mean stellar ages 
as well jWorthev & Ottaviani|[l997l) . Even though several studies 
have shown that the Balmer emission due to HII regions, AGN 
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Figure 5. The distribution of tlie (emission corrected) equivalent width 
(EW) for galaxies in our cluster sample, with line styles the same as in 
Fig.|4] As expected, the blue star-forming galaxies have the highest absorp- 
tion, showing the presence of < 1 Gyr old stellar population/s, while the 
red sequence galaxies show little evidence for absorption in H^. The blue 
passive galaxies have a broad distribution with equal probability of finding 
Hg in emission (negative EW) or absorption. The distribution for red star- 
forming galaxies is bimodal, which, in conjunction with the distribution 
of their internal extinction (Fig. |4), indicates the presence of two distinct 
galaxy populations. 




0^4000 

Figure 6. The distribution of the (Nelll continuum connected) D„4000 val- 
ues, with line styles the same as in Fig. |4] As expected, the red sequence 
galaxies have the largest values, while the blue star-forming galaxies the 
smallest values for Dn4000 respectively. The blue passive galaxies have 
a small offset (~ 0.2) from their star-forming counterparts, revealing the 
presence of slightly older stellar populations. The histogram for red star- 
forming galaxies has a statistically significant bimodality (see text), similar 
to that seen in Figs.|4]&ll] 

and/or planetary nebulae can fill in the underlying Ba lmer absorp- 
tion l i nes, causing age estim ates to be spuriously high jTrager et alj 
l200(]|;|Prochaska et alj|2007h . the EW remains a popular stellar 
age indicator because it is less affected than the lower o rder Balmer 
lines, due to the steep Balmer de crement in emission jOsterbro^ 
ll989l ; IWorthev & Ottavianill997h . 

Here we use the measured values of the (emission cor- 
rected) Hi EW, and the (Nelll continuum corrected) D„4000 from 



iKauffmann et alj ( l2003bl) . who, for t he latter, have a dopted a modi- 
fied version of the original definitio n teruzuall 19831) of the D„4000 
introduced by Balog h et alJ ( Il99 9^. A stronger D,i4000 is an in- 
dicator of a metal-rich inter-stellar medium, which, together with 
the Hf EW, indicates whether the galaxy has been forming stars 
continuously (intermediate/high D„4000 plus low H^ EW) or in 
bursts (low Dn4000 plu s high Rg EW) over the past 1-2 Gyr 
IKauffman n et al.ll2003bh . We note that the observational error in 
D„4000 index is very small, typically ~0.02, but the same is not 
true f or the Hg EW, which has a n uncertainty of 1-2 A in m easure- 
ment IKauffmann et al.ll2003bl) . IKauffmann et al J ( l2003bl) discuss 
the physical constraints leading to these uncertainties, and show 
that the H^ index can assume only a selected range of values for 
all possible star formation histories. The large observational errors 
on the Kg EW imply that the estimated value of this index is con- 
strained by the parameter space spanned by the acceptable models. 

The red sequence galaxies in our sample show H^ in emis- 
sion, and yield large values for the D„4000, confirming the fact 
that their main ingredient is a population of old, passively evolving 
stars, which determines the overall galaxy colours and dominates 
the spectra. The high observational uncertainty in the EW(Hj) im- 
plies that it is hard to determine whether H^ is actually occurring in 
emission, or the negative EWs are a result of a high correction fac- 
tor in individual galaxies. This is a common problem with the use of 
Hg index in individual systems, but this might not be an important 
issue for this statistical study. The blue star-forming galaxies, in 
sharp contrast, yield small values for the D„4000 and strong H^ in 
absorption. The blue passive galaxies follow the blue star-forming 
galaxies, with the mean value for the D„4000 systematically higher 
by ~ 0.2, but with a very broad distribution of H^ EW, revealing 
a mix of galaxies with a large scatter in the ages of their youngest 
stellar populations. In red star-forming galaxies, the distributions of 
Hg EW as well as the D„4000 are bimodal. The KMM test prefers 
a 2-mode fit over a 1-mode fit with likelihood ratios of 78 and 84, 
for Hg EW and D„4000 respectively, with zero probability for the 
null hypothesis in both cases. 

We also check whether this bimodality is a result of the thresh- 
old chosen to define whether a galaxy is star-forming or not. As for 
the internal extinction, the KMM test results remain statistically 
significant, for both Ha EW and the D„4000, when the red star- 
forming galaxies are selected with an even higher SSFR threshold 
(log SFR/M* = -10 yr~^), showing that these results are not sensi- 
tive to how the four categories are chosen, within reasonable limits. 
The above results become even more convincing when we plot the 
four different galaxy populations on the plane defined by the two 
indices D„4000 and EW(Ha) (Fig. |7ll. IKauffmann eTall ( l2003d) 
have shown that these two indices taken together can very effec- 
tively determine the star formation history of a galaxy. As shown in 
Fig- 13 this indeed seem to be the case. The blue star-forming and 
the passive, red sequence galaxies occupy the opposite ends of the 
space. But whilst the blue passive galaxies seem to follow the blue 
star-forming galaxies, although with somewhat lower EW(Ha) and 
higher D,i4000, the red star-forming galaxies span the entire range 
of values in both indices, as would be expected given their bimodal 
distribution. 

3.3 AGN and star formation 

Figs.|4l|5]|6|and|7]clearly show that while most of the cluster galax- 
ies follow the classical picture, blue colours correlating with higher 
SSFR, and redder colours with negligible star formation activity, 
there is a non-negligible fraction of galaxies in either class which 
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Figure 9. The central panel of this plot shows the {g — r)*^'^ colour of all the spectroscopically covered SDSS DR4 galaxies (grey dots) in the colour-SSFR 
space (same as Fig.|2j, found in the same redshift range and magnitude limit as our cluster galaxy sample. Over-plotted are 1.4 GHz radio sources from 
the NVSS (cyan), IRAS point sources (PSC, blue) and IRAS faint sources (FSC, red) respectively. The top and the right panel show the histograms of these 
quantities along the two axes. As expected, most of the radio-loud AGN concentrated along the red sequence, and IR bright sources in the blue star-forming 
galaxies. It is interesting to note that the red star-forming galaxies (top right quadrant, central panel), a large fraction of which are expected to be dust-reddened, 
do not seem to have many IRAS sources. This supports our conclusion here that a non-negligible fraction of the red star-forming galaxies obtain their colour 
from the presence of excessive metals. 



deviate from the norm. It is thus crucial to understand the proper- 
ties of these galaxies and identify the cause of their occurrence. In 
order to do so, we make use of the classification of galaxies, based 
on the BPT diagram, as star-forming, composite, AGN and those 
without emission lines (unclassified), as mentioned above in i]2.3l 

Ignoring, for the time being, the unclassified galaxies for 
which the nature of star formation activity cannot be conclusively 
determined from the spectroscopic data, we find that >30% of the 
red star-forming galaxies and -^50% of the red sequence galaxies 
are classified as AGN. This agrees with the well-known observa- 
tion that AGN are predominantly hosted by the massive elliptical 
galaxies. In Fig. |9] we plot 1.4GHz radio sources (from NVSS, 
iBest et al., | 2005 ) and IRAS 60^m sources found in the SDSS DR2 
I Cao et al.l2()06l) . along with all the galaxies found in the same red- 
shift range (z ^ 0.12) and magnitude limit (Mr ^ —20.5) in our 
galaxy sample. The radio-loud AGN are clearly seen to be hosted 
mostly in the red sequence galaxies. Interestingly, not many of the 



red star-forming galaxies are found in the IRAS 60^m source cata- 
logues, disagreeing with the usual expectation that the colour of the 
red star-forming galaxies is a result of the presence of dust around 
star-forming gas clouds. We further explore this issue in i]4.2| 

3.4 Aperture vs global colours 

An SDSS spectrum comes from a fixed aperture, which is equal 
to the diameter of the fibres (3") used in the spectrograph. Since 
this aperture would encompass a varying fraction of the light of a 
galaxy depending upon its size and distance, going up to ~40-60% 
at z = 0.1, the spectrum might not, in most cases, be a true rep- 
resentation of the galaxy's mean overall star formation activity. In 
this work, we compare photometric colours and spectroscopically 
derived specific star formation rates, which come from indepen- 
dent observations of the same galaxies. It would be interesting to 
see how the colour derived from the fibre aperture, as opposed to 
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Figure 7. The relation between the equivalent width of the Hs line EW(H|5) 
and the 4000A break index Dn4000, for all galaxies in our sample. These 
two indices together form a robust diagnostic of the star formation history 
of a ga laxy, irrespective of the stellar ages and metallicity iKauffmatm et alj 
l2003al also see text). The red sequence galaxies (grey scale) form a core 
in the passive region of the plane, while the blue star-forming galaxies 
(blue contours) occupy the region indicating recent and ongoing star for- 
mation. The blue passive galaxies (green contours) closely follow the blue 
stai' forming galaxies, but with somewhat higher Dn4000 and low EW(Hi ), 
indicating recent truncation of star formation. The innermost contours rep- 
resent 24 galaxies in both cases, and the number decreases by 3 at every 
successive level of contour The red star-forming galaxies seem to span the 
entire range in both indices defining this space, a relatively large fraction of 
them being passively evolving just like the red sequence galaxies, but a non- 
negligible fraction intruding the space dominated by the blue star-forming 
galaxies. 
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Figure 8. The cumulative distributions of g — r aperture colour, measured 
within the 3" SDSS fibre, relative to the overall colour of the entire galaxy. 
Atz~0.1, > 60% of the light in the fibre comes from the bulge and hence 
the fibre colour can be assumed to be a proxy for the colour of the bulge. 
The slope of the distributions for blue galaxies is much steeper than that 
of the red ones, showing the large spread in the relative colour, in'espective 
of the SSFR. Although the red galaxies have a relatively small spread in 
colour, a significant number of red star-forming galaxies have bluer bulges. 
This could be the signature of the presence of active nuclei supporting star- 
formation activity in the core of the galaxies, while the overall galaxy colour 
remains dominated by the older stellar population elsewhere in the galaxy. 



the value averaged over the entire galaxy, relates to SSFR. For in- 
stance, if an early-type galaxy has significant star formation in the 
core, we would expect to see a bluer core, while the overall photo- 
metric colour would be dominated by the older stellar populations 
elsewhere in the galaxy. 

In Fig [8] we plot the difference between the colour measured 
within the fibre and the overall colour of the galaxy, normalised by 
the overall colour, separately for galaxies belonging to the four cat- 
egories. A noticeable characteristic of this plot is that irrespective 
of their star formation properties, the blue and red galaxies have 
very different slopes for this cumulative distribution. The steeper 
slope for the blue galaxies is a consequence of their wider range in 
colour. While most of the galaxies have fibre colours that are red- 
der than the overall galaxy colours, ~15% of the blue star-forming 
and ~15% of the red star-forming and red sequence galaxies have 
bluer cores. This suggests the presence of centrally concentrated 
star formation in these galaxies. It must be noted that although rel- 
ative fibre colour is a good first approximation to the colour of the 
bulge of a galaxy, in this work where we stack galaxies over a range 
of redshift, the ratio of central and overall colour is not derived on 
the same physical scale for all galaxies. This result should therefore 
be considered to be a qualitative one. 



4 DISCUSSION 

Numerous studies of the evolution of galaxies have made use 
of either global photo metric proper ties such as colour (e.g. 
iButcher & Oemleill984l) , morphology ( lDressletll980l) or luminos- 
ity of a ga laxy, or spect roscopic properties such as circular ve- 
locity (e.g. iGraves et a l."2007). the D„4000, EW of H« absorp- 
tion line (e.g. iBa logh et a l.iil999l), emission lines of Hq or Oil 
(e.g. iLewis et al. ' 2002; G omez et al.l |2003| ; iBalogh et alj l2004bl: 
iHaines et alJl2007 ) or integrated SFR and/or SSFR (B04). Often 
samples are chosen based on colour or star formation properties. 
Are the properties derived from photometry and from spectra cor- 
related enough to yield similar results? 

Blue galaxies are thought to be of late-type morphology and 
and have significant star formation activity, while red galaxies 
are early-type and show little or no evidence of current star for- 
mation. However, with the recent availability of extensive multi- 
wavelength photometric databases, and of large optical spectro- 
scopic surveys, several authors (e.g. IWolf, Grav & Meisenheimetj 
'2005; Bildfell et al. 2008; Haines et al. 2008; Kovama et al. 200i 
i Saintonge^ Tran & Holden 2008; Graves, Faber, & Schiavon 20091; 
lwolfetalj|2009ir have suggested that optical colour may not nec- 
essarily be a good representation of a galaxy's evolutionary state. 
In this paper, we have explored both photometric colours and spec- 
troscopic specific star formation rate from optical observations of 
galaxies, in a wide range of environments, in and near rich clusters 
in the local Universe, to examine cases where the usual correlation 
between the two fails to hold. 



We find that >80% of the red galaxies ({g-r)°'^ > 0.8) that 
are classified to be part of the red sequence of clusters, on the basis 
of photometric colours, also have very little or no star formation 
(log SFR/M* 55-10.16 yr"\ see Fig.Hl. Similarly, 65% of the blue 
galaxies found within 3r2oo of rich clusters also have high SSFR. 
However, in addition to these 'normal' galaxies in which the colour 
and spectroscopic SFR/M* correspond well with each other, a sig- 
nificant fraction of galaxies deviate from the norm, and would be 
classified differently, if treated solely on the basis of one of their 
spectroscopic and photometric properties. 
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Figure 10. The distribution of the concentration parameter, C = Rgo/Rso 
for the red sequence (red solid) and blue passive (blue dashed) galaxies 
respectively. The thick and thin lines represent the distributions with and 
without unclassified galaxies in each case respectively. The blue passive 
galaxies clearly have late-type morphologies along with bluer colours. The 
fact that these galaxies have very low SSFR thus make them candidates for 
progenitors of the red sequence galaxies, in which star formation has been 
recently (< 1 Gyr) truncated, eventually leading to red colour and early- 
type spheroidal stiucture typical of the present day cluster red-sequence 
galaxies. 



4.1 Blue passive galaxies 

The blue passive galaxies in our cluster sample (Fig.|2} have spec- 
troscopic properties quite similar to their star-forming counterparts. 
The colour of a galaxy cannot be expected to, and indeed does not, 
correlate in the same way with environmental properties, as the 
star formation properties do (e.g., Kauffmann et al. 2003^. Our 
analysis shows that the values for the D„4000 (Fig. |6ll and in- 
ternal extinction (A^; Fig. |4j for blue passive galaxies are sim- 
ilar to those of blue star-forming galaxies, but the Hs EW dis- 
tribution (Fig. |5]( is very different for the two populations. Since 
our sample is chosen to be in or near rich clusters, the discrep- 
ancy between their mean photometric colour, and their SSFR 
derived from spectra, suggests two possible scenarios: (i) these 
are galaxies in which star formation has been recently (~l-2 
Gyr) shut off due t o the imp act of local and/or gl obal environ- 
ment (e.g.lPorter & Ravchaudhurvll2005l,l2007l:[ P OTter et al. 200a ; 



iMahaian. Ravchaudhurv & Pimbblelllin preparationh . or (ii) these 
are bulge dominated galaxies where the mean colour is dominated 
by the young stellar populations distributed in the disc but the spec- 
trum sourced from an almost dead centre shows little sign of star 
formation activity. 

A galaxy that has experienced a major burst of star forma- 
tion 1-2 Gyr ago will show significant absorption in H^, which is 
found in a non-negligible fraction of the blue passive galaxies here 
(Figs.[5]and|7](. This is because the A stars which significantly con- 
tribute to the Us absorption, have a lifetime of ~1.5 Gyr, but the 
burst of star formation produces enough ionizing radiation to fill 
in the Kg absorption line. But the absorption line shows up in the 
galaxy spectrum once the starburst has been truncated/ended. 

Fig.[TO]provides statistically significant evidence in favour of 
(ii). As a morphological index, here we use the concentration pa- 
rameter C = Rgo/Rso, where R^, is the Petrosian radius encom- 
passing x% of the galaxy's light in the r-band from SDSS images. 



Red & blue star-forming galaxies 

(including unclassified galaxies) 

Red & blue star-forming galaxies 

(excluding unclassified galaxies) 

Red star-forming galaxies 

(including & excluding unclassified galaxies) 

Blue star-forming galaxies 

(including & excluding unclassified galaxies) 

Red & blue passive galaxies 
(including unclassified galaxies) 
Red & blue passive galaxies 
(excluding unclassified galaxies) 
Red passive galaxies 

(including & excluding unclassified galaxies) 
Blue passive galaxies 

(including & excluding unclassified galaxies) 



0.67 0.0000 

0.60 0.0000 

0.08 0.1870 

0.02 0.9999 

0.35 0.0000 

0.29 0.0000 

0.08 0.1772 

0.01 0.9999 



In Fig. [To] we plot the distribution of C for the red and blue pas- 
sive galaxies respectively. The bl ue passive galaxies clearly show 
late-type morphologies (C < 2.6; IStrateva et al.ll200ll) . while their 
red counterparts mostly have values of concentration typical of 
spheroidals. The galaxies in the former category are similar to the 
passive blue spirals that have been studied by other authors (e.g. 
iGoto et al.ll20()3al ; Iwild et al.ll2008l) . 

By analysing the spatially resolved long-slit spe ctroscopy of a 
few pa ssive spiral galaxies selected from the SDSS. Ilshigaki et al] 
( |2007|) find that the Hs absorption lines are more prominent in the 
outer regions of the disc, while the nuclear regions show strongest 
values for D„4000. Their findings are consistent with the idea that 
the star formation in passive late-type galaxies ceased a few Gyr 
ago. Hence, the blue passive galaxies can be considered to be the 
progenitors of their red counterparts, in which star formation has 
been recently shut-off: these galaxies will eventually acquire red- 
der colours and early-type m orphologies, over a period that de- 
pends upon their environment jGoto et al.l2003d) . lKaimappan et al.l 
I2OO9) have recently identified, characterized and studied the evo- 
lution of galaxies which, though morphologically classified as 
E/SOs, lie on the blue sequence in the colour-stellar mass space, 
and suggest that these systems almost always have a bluer outer 
disk, and that they are a population in transition, potentially evolv- 
ing on to the red sequence, or have (re)formed a disk and are about 
to fall back in the class of late-types. 



4.2 Red star-forming galaxies 

Even more interesting are red galaxies with signs of active star for- 
mation in their fibre spectra. We find that 53.4% of the 539 galaxies 
in this class have significant emission lines on the basis of which 
their SSFR has been evaluated. These red star-forming galaxies 
are made up of at least two kinds of galaxies (chosen from sim- 
ple limits in colour and SSFR). One kind appears to be very simi- 
lar in photometric properties to the red sequence galaxies but with 
spectroscopically derived properties consistent with their late-type 
counterparts, while the other class has both photometric and spec- 
troscopic properties similar to the red sequence galaxies, yet has a 
high SSFR. 

The origin of the latter class is not clear It is likely that these 
galaxies are currently experiencing a starburst, which has started 
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very recently (<0.5 Gyr). Stellar populations resulting from such 
a starburst would be detected in spectral indices such as Ha EW 
(leading to high values of SFR/M*), but will not dominate the 
D„4000 and the H« EW for another ~0.5-l Gyr. The former sub- 
population is most likely to be bright cluster galax ies (BCGs) , 
with nuclear star-formation linked with the AGN (see lGotoll2006l : 
iBildfell et al.ll2008l : lReichard et al.ll2009l) . In the following we turn 
our attention to this sub-population. 



4.2.1 Star formation at the core of giant ellipticals 

The hypothesis that AGN and star formation activity in the core of 
a galaxy may be linked is supported by the fact that ~30% of the 
red star-forming galaxies with emission lines are classified as AGN 
on the BPT diagram. The fact that 50% of the emission-line red 
sequence galaxies (defined according to Fig. O are also classified 
as AGN does not contradict this scenario because either (i) the red 
sequence AGN galaxies in general have high stellar mass, and so 
can have low SSFR even though they have significantly high SFR, 
or (ii) they do not have enough cold gas to form stars, because of 
the presence of an active nuclei. 

A significant fraction of giant ellipticals show evidence of on- 
going star formation or signs of recent (<2 Gyr) star formation in 
their cores, particularly in environments where the li kelihood of re- 
cent mergers is high (e.g.„ McDermid et al. 2006; Nol an et alj2006l : 
iNolan. Ravchaudhurv. & Kaban 20070- For a s ample of BCGs in 48 



OTft.: 

X-ray luminous clusters. IBildfell et al.l ( |2008|) find that 25% of the 
BCGs have colour profiles (g—r) that turn bluer toward the centre. 
Our radial distribution of galaxies in this category (Fig. [3), com- 
pared to the distribution of blue star-forming galaxies, is consistent 

with this (see m.5\. 

Elsewhere lCjallazzi et ai]l l2009l) analyse an extensive dataset, 
covering UV to IR SEDs for galaxies in the Abell 901/902 cluster 
pair, to show that ~40% of the star-forming galaxies residing in in- 
termediate to high density environments have optically red colours. 
They show that these galaxies are not starbursts, and have SFRs 
similar to or lower than that of their blue counterparts. Whether the 
occurrence of blue cores in red sequence galaxies is linked to their 
nuclear activity, and whether the nuclear activity itself is modulated 
by the immediate environment of a galaxy (such as the presence of 
close neighbours), remains a debatable issue. We intend to return 
to this in future w ork inv olving a more detailed study of nuclear 
starbursts (also see iLi et a l. 2008; Re ichard et al.l2009h . 

The scenario discussed in (ii) above complements the blue 
passive galaxy population discussed in ^4.1\ - massive galaxies 
without much cold gas could be formed in dry major merg- 
ers of the passive blue progenitors, while high speed encoun- 
ters amongst them in moderately dense environments , such 
as in the infall regions of clusters (e.g. IPorter et al.l 1200 8l ; 
iMahaian. Ravchaudhurv & Pimbble^lin preparationl) . could lead to 
excessive loss of cold gas, leaving the galaxy to evolve passively 
thereafter. 

In a spatially resolved sp ectroscopic study of 4 8 early-type 
galaxies, the SAURON team jde Zeeuw et al.l l2002h simultane- 
ously observ ed the CO(l-O) and CO(2-l) lines with 28% detection 
rate (12/43) dCombes et al]|2007h . In this sample they find that the 
CO-rich galaxies show strongest evidence of star formation in the 
central (~4x4 arcsec^) regions and also tend to have higher and 
lower Fe5015 and Mgfo indices. The results presented in this work 
also support various studies based on the UV data obtained by the 
Galaxy Evolution Explorer (GALEX) which suggests that a signif- 
icant fraction (~15%) of bright (Mr < —22) early-type galaxies 
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Figure 11. The distribution of stellar metallicities (from iGallazzi et alj 
l2005i) for the red star-forming (red solid) and blue star-forming (blue 
dashed) cluster galaxies in our sample. The thick and thin curves are dis- 
tributions including and excluding the unclassified galaxies respectively in 
each case. The distributions evidently show that the stars in the red star- 
forming galaxies are much more metal-rich with respect to their bluer coun- 
terparts. Given that the internal extinction in galaxies belonging to the two 
categories does not show a significant difference (Fig.|4), presence of metal- 
rich stellai' populations in the red star-forming galaxies seems to be a hkely 
cause of their redder colours. 



show signature of recent (< I Gyr ol d ) star formation (e.gfYi et al.l 
l l2005h ; also see lCrocker etaP ( l2009h ; lGallazzi et alj l l2009ly K 



4.2.2 Stellar and gas phase metallicities 

Since star formation in a galaxy is often associated with the 
presence of dust, one would naively expect these red star-forming 
galaxies to be the dusty s tarburst galaxies widel y discussed in 
the lit erature JWolf, Gray & Meisenh eime^ |2005|; IBildfell et al.l 
2008 !; iKovamaetalJ 1200 8; Saint onge. Tran & Holdei] |2008|; 
IGallazzi et al.ll2009l ; IWolfet al..i2009l) . However, we find that only 
~ 50% of the galaxies in this class show presence of any optical 
extinction (Figs.|4]&|9]l. So why are the rest of the galaxies red? 

To address this question, we explore two other factors af- 
fecting galaxy colours, namely age and metallicity. We use stellar 
metallicities an d mean lumin osity-weighted stellar ages from the 
work of Gallaz zTet al.l ( l2005h . which are estimated from 150,000 
Monte Carlo simulations of a full range of physically plausible star 
formation histories. These models are logarithmically distributed 
in a metallicity range of 0.2-2.5 Zq. In a given model, all stars 
have the same metallicity, interpreted as the 'luminosity-weighted 
stellar metallicity'. For each model in the library, Gallazzi et ;d] 
( 2005) measure the strength of the D„4000, and the H/j, Hsa+UjA, 
[Mg2Fe] and [MgFe]' spectral indices in the same way as they are 
measured in the SDSS spectrum of a galaxy. By comparing the 
strengths of these indices in the observed spectrum to the model 
spectra in their library, the probability density functions (PDFs) of 
several physical parameters, such as age and metallicity of a galaxy, 
are constructed. 

We use the r-band luminosity weighted ages (hereafter re- 
ferred to as 'mean galaxy ages') and stellar metallicities derived 
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from these PDpfl For galaxies in our cluster sample, we find that 
the red star-forming galaxies are systematically older than their 
blue counterparts by 5-7 Gyr. In Fig.[TT] we compare the distribu- 
tions of stellar metallicities of the red and blue star-forming galax- 
ies. We confirm that the two distributions are statistically different 
by using the Kolmogorov- Smirnoff (K-S) statistics, which finds the 
probability of two datasets being drawn from the same parent pop- 
ulation. The K-S test finds the maximum difference between the 
cumulative frequency distributions of the red and blue star-forming 
galaxies to be D = 0.35, with the probability of the null hypothesis 
being satisfied to be essentially zero. 

Further support for the presence of excessive metals in the 
red star-forming galaxies is provided by a comparis on between the 

i ;as-ph ase metallicities (log[0/H]+12) estimated bv lXremonti et al.l 
2004Pl who have compiled a sample of ~53,000 star-forming 
galaxies from the entire SDSS DR4 (selected to have Hq,, & 
[Nil] emission lines at > 5cr detection). We find that the cumula- 
tive distributions of the gas-phase metallicities of the red and blue 
populations are statistically different. The K-S test yields i3 = 0.18 
at a significance level of 0.043. Though the two samples compared 
here have widely different numbers of galaxies, our conclusion is 
consistent with that of Ellison et al. (2008), who show that for a 
sample of 43,690 galaxies, selected from the SDSS DR4, the gas- 
phase metallicities of galaxies with high SSFR are lower by 0.2 
dex, compared to those with low SSFR at fixed stellar mass. The 
fact that their metallicities are estimated by a comparable but inde- 
pendent method (see Kewley & Dopita (2002} for details), further 
strengthens our result. 

Around 35% of all red star- forming galaxies are classified 
as star-forming or composite g alaxies from their spectra. This is 
in agreement with the results of [Graves. Faber. & SchiavonI ^2009} 
who also find ~30% of nearby red sequence galaxies from SDSS 
show presence of emission lines in their spectrum. They analyse 
LINER like emission in these galaxies to show that they are sys- 
tematically younger by 2.5-3 Gyrs than their quiescent counter- 
parts without emissio n at fixed rotational velocity. At z ~ 0.1, 
iTremonti et alj j2004h have shown that the metallicity is well corre- 
lated with the stellar mass of the galaxy. These quantities correlate 
well for 8.5^1og M* ^ 10.5 Mq, but the relation flattens thereafter. 
The fact that ^60% of our red galaxies have log M* in the range 
10.6-10.9 M0 supports the hypothesis that a significant fraction 
of the red star-forming galaxies obtain their red colours from the 
residual older generations of stars. 

We have also shown that ~50% of all the red star-forming 
galaxies are likely to be attenuated by >0.2 mag in the SDSS z- 
band (Fig.Q, and have bimodal Hf EW and D„4000 distributions 
(Figs.[5][6]and|7} which are very different from that of the typical 
red sequence galaxies. As discussed in Sj3T|and ^3.2\ the red star- 
forming class of galaxies remain bimodal in the distributions of Az , 
Us and D„4000, even if the selection criterion based on SFR/M* 
is changed to a significantly different value (e.g. ~15-percentile, 
corresponding to log SFR/M* = -10 yr^^). This clearly indicates 
that the component of red star-forming galaxies that are similar in 
properties to the red sequence galaxies are not merely a result of 
scatter from the red sequence (quadrant 1; Fig.|2ll. 



^ http://www.mpa-garching.mpg.de/SDSS/DR4/Data/stellarmet.htmI 
^ http://www.mpa-garching.mpg.de/SDSS/DR4/Data/oh_catalogue.html 



4.3 Implications for the Butcher-Oemler effect 

iButcher & Oemleil j 19841) found that clusters at moderate to 
high redshift contain an excess of blue galaxies, compared to 
their low-redshift counterparts. For the redshift regime of our 
sample (z ~ 0.1), they found a uniform blue fraction f(, of 
~0.03 in all clusters within R30, the radius containing 30% of 
the cluster's red sequence population. This exercise has been 
repeated several times in recen t years for many samples of 
clusters (|Butcher & Oemler 1984: Margon iner & Carvalhol I2OOOI; 
lEUingson et al. 2001 ; De Propris et al. 2004). Most studies of this 
kind define 'blue' galaxies in terms of a broadband colour that is 
bluer by 0.2 mag than that of the cluster's red sequence galaxies. 
From the above discussion, it is clear that this cut leaves out a sig- 
nificant fraction of star-forming galaxies that have redder broad- 
band colours, and classifies some passively evolving galaxies as 
'star-forming'. 

iMargoniner & Carvalhol ( I2OOOI) found the blue fraction of 
galaxies to be fi,=0.03±0.09, similar to lSutcher & OemleJ ( Il984l) . 
for 44 Abell clusters at 0.03^z5C0.38 within a fixed cluster-centric 
distance of 0.7 Mpc, which translates int o the mean R30 for mos t 
Abell clusters below z < 0.1. However. IMargoniner et al.l ( l200lh 
measured the fraction to be (1.24 ± 0.07)2 — 0.01 for clusters at 
z ^ 0.25, again for galaxi es within a fixed cluster-centric aper- 
ture o f 0.7 Mpc. Elsewhere dEUingson et al.ll200ll : lDe Propris et al.l 
|2004 . various multiples of cluster-centric radius, scaled with r2oo , 
have been us ed to show the effect o f chosen aperture size on the 
blue fraction. lEllingson et al.l ( l200lh suggest that the origin of the 
Butcher-Oemler effect lies in the fact that the relative fraction of 
'blue' galaxies on the outskirts of the clusters at higher redshifts 
is higher than in their local counterparts. If confirmed, this result 
can have a critical impact on the study of evolution of galaxy prop- 
erties. It is worth mentioning here that even though most of these 
studies use optical spectra for most or all of their sample (except for 
iButcher & OemleJ ( Il984h ). the spectroscopic information is used 
only for assigning cluster membership. 

Recently, with the increasing availability of large multi- 
wavelength datasets, the possibility of the Butcher-Oemler effect 
being observed in other fundamental galaxy properties (e.g. mor- 
phology. Goto et al. 2003b), or in non-optical data (e.g. mid-IR, 
Saint onge. Tran & Holden 2008), is being explored. The study 
of iWolf. Grav & Meisenheimer. C2005.) unraveled the presence of 
young stars and dust in the red seque nce galaxies in a cluster pair 
at z~0.17, while lBildfell et al.l ( l2008h find that the colour profile of 
25% of the brightest cluster galaxies (BCGs) turn bluer by 0.5-1.0 
mag towards their centres, compared to the average colours seen in 
the cluster's red sequence. All of this has serious implications for 
Butcher-Oemler like effects based on global photometric proper- 
ties. 

Our sample provides some additional insight to this discus- 
sion. The presence of metals, and of colour gradients related to 
SFR in a galaxy, will affect the estimated fraction of 'star-forming' 
galaxies in low redshift clusters. The fact that nearby cluster galax- 
ies are more metal-rich than those in high-redshift clusters, will 
artificially enhance the fraction of passive galaxies that are selected 
solely on the basis of their colour. Such an effect is usually not 
accounted for while comparing the blue/star-forming galaxy frac- 
tion in clusters at different redshifts. A quantitative analysis of this 
effect requires a consistent study of an unbiased sample of mutu- 
ally comparable galaxy clusters spanning a wide redshift range, for 
which both photometric and spectroscopic multi-wavelength data 
is available. We leave this for a later consideration. 
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5 SUMMARY AND CONCLUSION 

We divide a sample of >6,000 galaxies found in or near rich Abell 
clusters (z^0.12) in the SDSS spectroscopic catalogue, into four 
populations, defined by placing simple limits in {g — r)^'^ colour 
(derived from photometry) and specific star formation rate (derived 
from the detailed modelling of optical spectra). While most blue 
galaxies have evidence of star formation, and red galaxies do not, 
this identifies two significant populations of blue passive galaxies 
and red star-forming galaxies. 

We trace the spectroscopic and photometric properties of 
galaxies in these four sub-samples in detail, using data available 
from SDSS DR4. The main results of our analysis can be sum- 
marised as follows: 

• For the blue passive galaxies, the dust content and values for 
the D„4000 are similar to those of the blue star-forming galaxies, 
indicating similar metal content, mean ages and internal extinction. 

• Since the blue passive galaxies have a broad distribution of Ha 
EW, along with late-type morphologies, they appear to be similar to 
the passive blue spirals studied by other authors. Our results support 
the hypothesis that these galaxies could be the progenitors of the 
red sequence galaxies in which star formation has been recently 
shut off due to the impact of local/global environment. 

• An intriguing find of this study is that only ~50% of the red 
star-forming galaxies have significant values of extinction, indica- 
tive of the presence of dust. 

• By comparing the stellar metallicities, gas-phase metallicities 
and mean ages of the red and blue star-forming galaxies, we find 
statistically significant evidence in support of the fact that, for at 
least half of the red star-forming galaxies, the red colours result 
from a relatively metal-rich stellar population. Our argument re- 
mains valid even if the criterion for classification as star-forming 
based on SSFR is changed within reasonable limits (15-50 per- 
centile, instead of the 25-percentile used in this work). 

• A significant (~15%) fraction of the red star-forming galaxies 
have blue cores, indicating recent star formation which might be 
linked to the AGN activity in these galaxies. 

• The population of passive red sequence galaxies increase to- 
ward the cluster centre, closely followed by the red star-forming 
galaxies, except that the latter does not rise steeply in the cluster 
core. 

• The blue star-forming galaxies are distributed uniformly at 1- 
3r2oo from the cluster centre, but decline steeply within r2oo • The 
blue passive galaxies on the other hand show a constant distribu- 
tion, except for a small increase around 2r2oo ■ 

• Galaxy evolution studies which base their sample selection 
criterion only on broadband colours to test Butcher-Oemler like ef- 
fects, may be significantly affected by the presence of 'complex' 
galaxy populations such as the red star-forming and blue passive 
galaxies discussed here. 
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